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In condensed-matter, level statistics has long been used to characterize the phases of a disordered 
system. We provide evidence within the context of a simple model that in a disordered large-N 
gauge theory with a gravity dual, there exist phases where the nearest neighbor spacing distribution 
of the unfolded pole spectra of generic two-point correlators is Poisson. This closely resembles the 
localized phase of the Anderson Hamiltonian. We perform two tests on our statistical hypothesis. 
One is based on a statistic defined in the context of Random Matrix Theory, the so-called A3, or 
spectral rigidity, proposed by Dyson and Mehta. The second is a ^-squared test. In our model, the 
results of both tests are consistent with the hypothesis that the pole spectra of two-point functions 
can be at least in two distinct phases; first a regular sequence and second a completely uncorrelated 
sequence with a Poisson nearest neighbor spacing distribution. 
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Real-world condensed-matter systems are always im- 
pure; realistic samples arc often described by non- 
translationally invariant interactions, about which only 
limited statistical information is available. As usual, 
perturbation theory has a limited range of applicabil- 
ity and answering questions concerning the phase struc- 
ture of a disordered system generally requires knowledge 
of non-perturbative dynamics. A prime example is the 
celebrated Anderson Localization (AL) phenomenon 
which, in sufficiently high dimensions, is entirely inacces- 
sible by perturbation theory. 

In recent years, the Anti-de Sitter/Conformal Field 
Theory (AdS/CFT) correspondence [2j or more gener- 
ally gauge/gravity duality has taught us a great deal 
about certain strongly correlated large-N gauge theo- 
ries 1 , by offering a dual gravitational perspective. It has 
provided us with a computational edge, in uncomfort- 
able regions of parameter space of field theories, where 
the conventional pcrturbativc techniques fail. Over the 
past few years, understanding quantum critical points 
in various condensed matter systems using gauge-gravity 
duality Q, 01, nicknamed AdS/Condensed-Matter The- 
ory (AdS/CMT), has emerged as a new avenue of re- 
search, where the correspondence, at least qualitatively, 
has proven insightful. In this approach one writes down 
phenomenological gravitational dynamics which are ex- 
pected to capture some aspects of the condensed matter 
system under study. Due to technical challenges, disor- 
dered systems have not received much attention. In Q, 



the effects of perturbative disorder on transport phenom- 
ena were first pioneered. Authors of Q, 0| studied renor- 
malization group properties of disorder in holographic 
settings. Holography and structural disorder/glassiness 
were studied in Q, @- See also [To[ for a replica-trick 
formulation of disorder-averages in AdS/CFT. 

In this paper we seek to address the following press- 
ing question: what constitutes a "generic" hallmark of a 
disorder-induced phase transition in a system possessing 
a gravity dual? 2 A disorder-induced phase transition of 
central importance is the AL transition pj. In dealing 
with disordered systems in condensed-matter, one ap- 
proach has been to study statistical properties of the lev- 
els, i.e, the discrete energy spectrum of a given disordered 
system. The particular statistic of interest is the so-called 
nearest neighbor spacing distribution (NNSD), defined 
in the context of Random Matrix Theory (RMT) . Level 
spacing distribution is a reliable order-parameter to dis- 
tinguish localized from delocalized phases, when a disor- 
dered system undergoes an AL transition , [ljj , [H[ . 
For instance, the disorder-driven Metal-Insulator Tran- 
sition, induced by AL, is marked by a non-perturbative 
change in the NNSD of levels from Gaussian Orthogonal 
Ensemble (GOE) universality class 3 in the delocalized 
regime to Poisson in the localized phase [l4[ . 

In this paper we demonstrate, in the framework of a 



1 Here "N" denotes the number of colors. In this work, we 
working only in the gravity limit, where N is large. 
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2 It would be advantageous if the existence of the order-parameter 
does not rely on perturbation theory. 

3 Depending on microscopic symmetries of the Hamiltonian, the 
GOE here must be replaced by one of the other Gaussian En- 
sembles. 
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simple model, the existence of a similar phenomenon, 
where the "condensed matter system" is a field theory 
with a gravity dual and non-translationally invariant in- 
teractions. We work with quenched disorder only. 

At large-N, even at finite volume, level statistics is not 
a good and/or convenient observable. In this paper, we 
put forward an alternative proposal, suggesting that for 
a large-N gauge theory a candidate order-parameter is 
the NNSD of the pole spectra of generic two-point cor- 
relators. It is important to note that the discrete pole 
spectrum does not represent the full state content of a 
many-body system. We consider the case where the pole 
spectrum consists of at least one discrete pole sequence, 
containing large number of poles. In the toy model dis- 
cussed in the present work, there are infinitely many dis- 
crete poles all lying on the real axis. 

In a generic disordered system, there exists a set of 
operators with position-dependent source terms, 

deforming the theory. Whenever a gravity dual is avail- 
able, the operators {0*,} will be dual to bulk fluctua- 
tions with prescribed position-dependent boundary 
conditions. In general, the back-reacted bulk geometry 
will have no Killing symmetries. Using AdS/CFT, we 
compute the real-time two-point correlator of a generic 
composite operator 0\ in the boundary theory using this 
deformed bulk spacetime [l^]. The correlator is defined 
as 



(1) 

where A is the bulk field dual to 0\. Note that the two- 
point function poles encode the spectrum of states (in- 
cluding bound states and resonances) which could be cre- 
ated/annihilated by 0\. From this view point, it paral- 
lels the energy levels in the single-body Anderson Hamil- 
tonian problem of > where the use of level statistics 
was advocated. 

Our static bulk geometry is asymptotically global 
AdS 3 
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where p £ [0, §], ip e [0,27r]. The AdS radius is set 
to one. The AdS3 asymptotic boundary condition is 
implemented by requiring ^,^ p ^s. ~ 
thus allowing for non-trivial boundary metrics. The La- 
grangian density for the bulk field A is 
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where A is taken to be dual to a relevant deformation. 
We will only need the quadratic part of ([3]). Using the 
AdS/CFT dictionary, this setup is expected to capture 
aspects of a disordered large-N gauge theory defined on a 
1+1-dimensional space with general metric. The spatial 
direction has a circular topology. We call the limit where 
$ = * = "clean" . 



In order to remain general, we avoid limiting ourselves 
to any particular AdS / CMT embedding of the disordered 
boundary theory. Our attitude will be that although 
the background © solves AdS3 gravity coupled to mat- 
ter equations of motion, nature of the matter source is 
left unspecified. Nevertheless, we demand that the bulk 
spacetime satisfies the Null Energy Condition (NEC) and 
is free of singularities. 4 This is done by checking if 
Ra>p(, a (,^ > everywhere in the bulk for an arbitrary 
null vector ( a [16j]. This space of geometries could be 
called potentially realizable backgrounds (PRB). The im- 
portant task of complete characterization of the PRB is 
left for future explorations. To begin, we will be content 
to study a class of backgrounds, where only $ is turned 
on. For this class, equation of motion for the probe field 
A is the following partial differential equation (PDE) 

sinV>+d>+[l+^~ A t We' 4 * sin 2 p-/C($)]^ = 0, 



4cos 2 p 

where A = tan _1 ' 2 p e~' s 'ipe~ lulT and 



(4) 



= sm 2 p[d 2 p $ + (d p $) 2 ] + tanp 8 P <S> (5) 

+[%$ + (d v *) 2 ]. 

The problem of extracting poles of the two-point corre- 
lator ([T]) is the same as solving Q subject to boundary 
conditions as a PDE eigenvalue problem. To tackle the 
singular PDE in Q , only numerical techniques are avail- 
able. In the simulations, varying m\ did not cause any 
qualitative change in the results reported here. Through- 
out this paper, we set to 2 = —0.75. With this choice of 
mass, there are two consistent quantizations of the bulk. 
We pick the one in which conformal dimension A\ of the 
operator 0\ is 3/2. The imposed boundary conditions 
are regularity of A at the origin of the global coordinates 
and vanishing of ip at the boundary. A family of back- 
grounds ^ 
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l + /(p)sinV 9 = 0, 



(6) 



where f(p) is regular at p = ~ and f(p) ~ p 2 near 
p = to ensure the absence of singularities, contains 
spaces satisfying the NEC. It is necessary that f(p) has 
no extremum. 

Before proceeding, let us briefly define the NNSD in 
the context of RMT 17( . To meaningfully compare the 
fluctuation properties of two sequences of data, or dif- 
ferent segments of a given sequence 5 , it is necessary to 
first unfold the sequence. This ensures that the non- 
universal differences, reflected in a non-constant average 
level density are removed. The unfolding procedure sets 



4 Once temperature is introduced, a horizon boundary condition 
will replace this somewhat strong condition of no singularity. 

5 This could be any sequence: energy levels of complex nuclei or 
spectrum of the operator in Q, for instance. 
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the average level density to one. In practice, to unfold 
a given spectrum one needs to numerically fit the level 
density histogram and then Xi = J e * £(e)de, where {x{\ 
form the unfolded sequence, marks the beginning of 
the sequence, and £ denotes the fitted local level density. 

As stated before, the order-parameter of interest to us 
is the NNSD of the eigenvalue sequence of (j4]), defined 
in the following. Take {ei, - • • ,e m } to be a sequence 
of m numbers, say eigenvalues of a large random ma- 
trix. Denote the corresponding unfolded sequence de- 
fined above, as {xx, ■ • • ,x m }, where the sequence has 
further been sorted in an ascending order. Now form 
Si = (xi+i — Xi)/ (s), where the average is over the entire 
sequence {sj}. The histogram associated to the set {s.j} 
is the NNSD. Before delving into the numerics, a com- 
ment is in order. The NNSD is usually calculated for a 
data set collected from many realizations of the quenched 
disorder. Within the approach here, it is not clear that 
averaging over an ensemble of bulk geometries is the same 
as disorder-average in the boundary theory. We only fo- 
cus on individual realizations. We hope to come back to 
this in the future. 

Powerful spectral techniques are utilized to numeri- 
cally compute the spectrum of accurately [T^j. We 
Chebyshev expand the wave function tp in the radial 
direction. Along the spatial boundary direction, a 
Fourier expansion is performed. To completely discrctize 
the PDE, we perform collocation at Chebyshev-Gauss- 
Lobatto points In the "clean" limit, poles of the 

two-point function ([1]) are given by uj = ±(A>, + 2k + £), 
where k, t = 0. 1. • • ■ We use eigenvalues of the 

clean system as a check on the performance of our code. 
We refer to the collection of lowest eigenvalues com- 
puted numerically, and sorted in the ascending manner as 
our "data set/sequence". It is empirically observed that 
about one half of the total number of eigenvalues, pro- 
duced in a computer run, arc fully resolved. Furthermore, 
to ensure the accuracy of the computed eigenvalues, all 
the analysis below is repeated with higher spectral resolu- 
tion until convergence is observed. In practice, roundoff 
errors are the limiting factor for large simulations with 
high spectral resolution. The eigenvalues come in pairs, 
due to symmetry r — > — r. Without loss of generality, 
the positive branch is picked. 

Now we are prepared to present our numerical re- 
sults. In order to illustrate the central idea proposed 
in this paper, we study a simple model of the form (|6|), 
parametrized as f(p) = Wp 2 (l + p 2 ), satisfying the NEC. 
The parameter W > will control the "strength" of the 
deformation Note that W can not be gauged-away. 
The use of Statistics is justified due to large number of 
poles involved in the analysis. We give sample values of 
W which correspond to Poisson and non-Poisson phases 
of the NNSD of the pole spectra. A more systematic 
search on the space of PRB will be reported elsewhere. A 
common method of testing a statistical hypothesis about 
a given data set is to apply a ^-squared test. For a given 
value of W , we perform the x~ s °L uare< i test to compare 



the numerically computed NNSD of the eigenvalues of (|4]) 
against Poisson distribution. The result is summarized 
in the so called p-value index. The index determines 
whether there is statistical evidence against the null hy- 
pothesis. 6 If the associated p-value of the x-squared test 
is larger than a (commonly accepted) threshold value of 
0.05, data is considered consistent with the distribution. 
The level density plot for W = 1.44 appears in Fig. Q] 
There is no universal recipe for choosing the "right" num- 
ber of bins, although there exist recommended values for 
the bin numbers, if one has prior knowledge about the 
distribution. We avoid any bias in choosing the number 
of bins. A cautionary note is that one should not at- 
tempt to fit the histogram too well. In particular, a fit 
involving high degree polynomials is undesirable. We ap- 
ply linear regression to the level density histogram. The 
fitting line is also overlaid on the level density histogram 
in Fig. [TJ The corresponding histogram for the unfolded 




FIG. 1: (color online). Level density plot for W = 1.44. Total 
number of participating eigenvalues is Nd — 406. 



data is given in Fig. [3J The NNSD for W — 1.44 is pre- 
sented in Fig. [2] The corresponding p-value is computed 
to be approximately 0.93, which indicates that our data 
set is consistent with the null hypothesis. The NNSD 
histogram for W = 0.009 is shown in Fig. @] with a p- 
value «0, which indicates a non-Poisson NNSD. In the 
following, we turn to W = 1.44 case and perform yet an- 
other statistical tests to build confidence that the NNSD 
is indeed Poisson. In the context of RMT, Dyson and 
Mehta suggested several statistics for studying the spec- 
trum of random matrices, the more useful one being A3, 
the so-called spectral rigidity [2(| ■ This statistic is sensi- 
tive to the long-range correlations in the spectrum. It is 



6 The null hypothesis here is "for a given value of W , the corre- 
sponding NNSD is Poisson". 
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FIG. 2: (color online). The NNSD for 406 eigenvalues of 
the system (|4]) with W=1.44. The cyan curve is the Wigner 
surmise for the GOE. The red curve is Poisson distribution 
P(s) = exp(— s). 



FIG. 4: (color online). The NNSD for N d = 406 eigenval- 
ues of the system Q with W = 0.009. The cyan curve is 
the Wigner surmise for the GOE. The red curve represents 
Poisson distribution P(s) — exp(— s). 
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FIG. 3: (color online). Level density of the unfolded spectrum 
for W = 1.44. It is nearly a constant (normalized to one) with 
the expected statistical fluctuations. 



defined as follows 



A 3 (L) = z (min J [ 



[M{x) - Ax - B] 2 dx) = (Aj(L)>, 

(7) 

where the average is over V , the position of the begin- 
ning of the sequence, and M(x) is the integrated level 
density for the unfolded sequence. The symbol "min" 



indicates that As(L) is calculated for the values A and 
B for which the right hand side in minimized. It is well 
known that for a Poisson distribution of levels, for large 
values of L, the spectral rigidity has a linear behavior 
with a constant slope 



1 



m - —l. 



(8) 



Applying this to our data set of eigenvalues, the plot 
in Fig. [5] is obtained. As is clearly observed, the case 
W = 1.44 corresponds to a Poisson sequence, while 
W = 0.009 is non-Poisson. Within the approach taken 
here, it might not be meaningful to determine the criti- 
cal value W c where the transition happens, since it is not 
immediately clear how to connect W to disorder in the 
boundary field theory. All that is known is that W = is 
the "clean" limit and W ^ should correspond to "dis- 
ordered" . The observed saturation of the linear behavior 
is well understood from the work of M. V. Berry 21 1. 
We also applied our analysis to the case where the back- 
ground (J2J) has a more general (^-dependence. Moreover, 
we studied (unsystematically) the case where f(p) —> 
as the boundary was approached. For these cases, the 
results were qualitatively similar to the ones reported in 
this paper. 
Discussion- 

We have demonstrated that at least sub-sectors of the 
spectrum of a disordered large-N gauge theory with grav- 
ity dual can become entirely uncorrelated/random. It 
would be interesting to further characterize the transi- 
tion. It is also crucial to investigate how generic our 
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FIG. 5: (color online). The blue and red diamonds represent 
As(L) computed from the data set of Nd = 406 eigenvalues 
for W = 1.44 and W = 0.009, respectively. The line y = L/15 
indicates the expected behavior of the statistic if the sequence 
is Poisson distributed. 



the future is to find out whether change in the NNDS 
is correlated with other signatures of the AL transition 
observed in transport, for instance. 

The issue of dependence on the realization needs fur- 
ther clarification. Perhaps the right approach is to av- 
erage over many spectra, each generated by a particu- 
lar realization of the disordered system. It would be in- 
teresting to generalize our results to higher-dimensional 
settings. Simulations involving a horizon in the bulk will 
teach us about statistics of resonances in the thermal 
plasma of a disordered gauge theory. The corresponding 
results will be reported elsewhere. 

Another future direction is to specialize to a particular 
embedding of the disordered boundary theory, where one 
solves for the back-reacted gravity background and com- 
putes the critical disorder strength at which transition in 
the NNSD takes place. 
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